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We adapt the notion of processes to introduce an abstract framework 
for dynamics in finite time, i.e. on compact time sets. For linear finite- 
time processes a notion of hyperbolicity namely exponential monotonicity 
dichotomy (EMD) is introduced, thereby generalizing and unifying several 
existing approaches. We present a spectral theory for linear processes in 
a coherent way, based only on a logarithmic difference quotient. In this 
abstract setting we introduce a new topology, prove robustness of EMD 
and provide exact perturbation bounds. We suggest a new, intrinsic ap- 
proach for the investigation of linearizations of finite-time processes, in- 
cluding finite-time analogues of the local (un)stable manifold theorem and 
theorem of linearized asymptotic stability. As an application, we discuss 
our results for ordinary differential equations on a compact time-interval. 

1 Introduction 
1.1 Motivation 

In this article we deal with what is informally referred to as finite-time dynamics and 
is usually represented by ordinary differential equations on compact time-intervals, i.e. 

x = f(t,x), 

where, for instance, / 6 C(I x 1", R") and I = [t_, t+], t_ < t + . The need to analyze 
such equations arises in many applications such as transport problems in fluid, ocean or 
atmosphere dynamics (see for instance [2 7) for a recent review), but also increasingly 
in biological applications (see [H US]). There are at least two reasons why one is 
interested in dynamics on bounded time-sets: one is the interest in transient behavior 
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of solutions although the differential equation might be given on the whole real (half- 
)line, and the other one is the simple fact that the right hand side / is given only on 
a bounded time-set, e.g. when it is deduced from observations or measurements. In 
any case, classical, i.e. asymptotic, concepts do not apply directly to the finite-time 
situation. During the last years efforts were made to establish finite-time analogues to 
asymptotic notions such as hyperbolicity of trajectories [23] and linearizations [7ll28j. 
Lyapunov exponents [22 [311 123] and stable and unstable manifolds [231 [211 EE [5] . 

The present article contributes to the investigation of hyperbolicity and its implica- 
tions in finite time, develops the theory in a coherent way and is organized as follows: 
after introducing notations in [subsection 1.2 we propose an abstract framework for 
dynamics on compact time sets I, not necessarily intervals, in |section "2[ For that 
purpose, there is no need to start the investigation with solution operators generated 
by ordinary differential equations (ODEs), which are given on an interval containing 
I. Instead, our framework is based on the well-known notion of processes, thus we 
state explicitly all required regularity conditions. To the best of our knowledge, such 
a first-principles approach to finite-time dynamics is proposed for the first time. It 
prepares the ground for numerical analysis. Starting from a slight modification of the 
elementary notion of a difference quotient, namely the logarithmic difference quotient, 
we introduce growth rates as supremum and infimum over the logarithmic difference 
quotient of trajectories. The step-by-step procedure enables us to prove new continuity 
results, see Proposition 2.12 and [Proposition 2.16 Another new feature is the intro- 
duction of a natural topology on linear finite-time processes, which is consistent with 
the supremum metric on continuous linear right hand sides of ODEs. The topology 
allows to consider perturbations of linear processes as in |Corollary 2.17[ |Lemma 2.21| 
and Proposition 3.7 In section 3 we develop a spectral theory for linear finite-time pro- 
cesses which is essentially known for solution operators on time-sets consisting of two 
points [28j and for intervals [Jj 114] . The underlying hyperbolicity notion is motivated 
by several approaches [23 [7J [2H E3] ■ These are generalized and unified altogether. 
We call this kind of finite-time hyperbolicity suggestively exponential monotonicity 
dichotomy (EMD). In subsection 3.3 we make use of the introduced topology to estab- 
lish the robustness of EMD easily. In section 4 we introduce formally linearizations of 
finite-time processes along trajectories and establish finite-time analogues of classical 
linearization theorems. We integrate the notion of stable and unstable cones firstly 
introduced in [14] . but characterize them by bounds on the growth rate functions, and 
give a direct and intrinsic proof of a Local Stable and Unstable Manifold Theorem 
as well as a finite-time analogue of the Theorem of Linearized Asymptotic Stability. 
In subsection 4.3 we investigate the local relationship between a process and its lin- 
earization, which can be interpreted as a finite-time Hartman-Grobman-like Theorem. 
In |scction 5] we apply the general results of this article to ordinary differential equa- 
tions on compact time-intervals. It turns out that the concept of finite-time Lyapunov 
exponents (FTLEs) coincides with the EMD-spectrum on time-sets consisting of two 
points. Consequently, our definition of spectrum can be used for a generalization of 
the FTLE-concept to arbitrary compact time-sets. 
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1.2 Preliminaries and Notation 

Throughout this article IcE denotes a compact subset of the real numbers, where 
t_ := mini and t + := max I. We will have occasions to mention the set of all ordered 
pairs of numbers in I with unequal components. In accordance with the notion of a 
relation, we denote by 7^nxi : = {(t, s) £ I x I; ( ^ s} c I x I the unequal-relation on 

1 x I. 

For a set M we write idM for the identity function on M and 2 M for its power set. 
To distinguish between the evaluation of a function f: A — > B a,t x E A from the 
image of a subset A' C A we write /(x) for the former and f[A'} for the latter. 

In the following, we consider dynamics on the Banach space (R™, |-|) with an arbi- 
trary Banach space (vector) norm |-|, which in turn induces an operator norm denoted 
by ||-||. By S we denote the unit sphere with respect to the given Banach space norm 
on R". At some point we will make use of the natural Euclidean structure of R™ to 
refer to orthogonality. Nevertheless, we do not need specific Hilbert space arguments. 
As usual, we denote by L(R n ) the set of linear operators on R™ and by GL(n, R) the 
subset of invertible operators on R n . 

For a Lipschitz continuous function / we denote by |/| Lip the Lipschitz constant 
of /. For time-dependent functions, whether vector- or operator-valued, the notation 
IHloo denotes the supremum norm over time with respect to the respective norms. For 
a continuously diffcrcntiable function / : Bq x . . . x Bk — > B, where B, Bq, . . . , B). are 
(open subsets of) Banach spaces, we denote the partial derivative of / with respect 
to the j'-th variable evaluated at x e B n x . . . x _B fc by djf(x). Note that the first 
argument has index 0. In case / has only one argument, the derivative is also denoted 
by f'(x) or f(t) if the argument is the time variable. 

We write B(x,d) and B[x,S], respectively, for the open and closed ball around x 
with radius 5 € R>o- 

2 Finite-Time Processes & Growth Rates 

2.1 Processes 

The notion of a process or two-parameter semi-flow was originally introduced by Dafer- 
mos however, with slightly different conditions than we require. In particular, we 
include invertibility already in the definition. 

Definition 2.1 (Process, linear process, smooth process). We call a continuous func- 
tion 

<ys: I x I x R™ — > R™, (t,s,x) i-> (p(t,s,x), 

an (invertible) process on I with state space R n if it is Lipschitz continuous in the 
first argument and if for any t,s,r S I and x € R™ we have ip(t,t,-) — id^n and 
ip(t, s, •) o ip(s, r, ■) = (p(t, r, •). We denote by V(l, R n ) the set of invertible processes on 
I with state space R™. We call an invertible process $ € "P(I, R n ) a linear (invertible) 
process on I if for any t, s € I we have $(i, s, •) € GL(n, R) and the function t i-> 
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&(t, s, •) <E L(R™) is Lipschitz continuous with respect to the operator norm. To 
emphasize the linearity in the last argument we will write <I>(i, s)x instead of s, x) 
for t, s € I, x € R n . We denote by £7 7 (I, R n ) the set of linear invertible processes on I 
with state space R™. Let k 6 IN>o. We call an invertible process ip a C k -process on I if 
for any t, s e I and x e R" we have y>(i, s, •) € C fc (R™, R") and d 2 cp(-,s, x) € ^(R™) 1 
is Lipschitz continuous. 

In the following, we always consider invertible processes and hence skip the word 
invertible. The required Lipschitz continuity of t i— > $(£,<_) in Definition 2.1 for a 
linear process $ is rather a technical assumption than an integral part of the notion 
of a process. In principle, it could be replaced by the weaker, but again technical, 
assumption of absolute continuity. For convenience, we will assume Lipschitz continu- 
ity in the sequel and point out which implications the weaker assumption would have. 
However, the necessity for a stronger continuity assumption than just continuity will 
become apparent in the course of this article. 

We start our investigations with linear processes on I. 

Lemma 2.2. Let $ E CV(I,R n ). Then the function 

IxI^L(E n ), (t,s) h-> $(t,s), 

is uniformly continuous with respect to the (induced) operator norm. 

Proof. This is a direct consequence of the uniform continuity of ^| IXIxB rQ jr D 

By the compactness of I we obtain in the next lemma the continuous dependence of 
the norm of trajectories on the initial value, uniformly in time. 

Lemma 2.3. Let ip £ P(I,R") and $ e CV(I,R n ). Then the function 

f: R" -> C(I,R"), i^(i4^,L)i), 

is continuous and the functions 

g: R" — s- C(I, R"), x ^ (t h-> t-)x) , 

h: R" -> C(I,R >0 ), ojh. (t h^- |$(t,i_)a;|) 

are Lipschitz continuous. 

Proof. To see the continuity of / in some x € R n it suffices to restrict the continuous 
function ip to the compact set I x I x B[x, 1] C R 2+ ™, thereby turning it into a 
uniformly continuous function. The continuity of / in x follows directly from that 
uniform continuity. By proving Lipschitz continuity for h we prove the assertion for g 
along the way. Let x,y S R™ and estimate 

IIM*) - Hv)\L - - M-^yiw^ 

<sup{|$(t,t_)(a:-y)|; t e 1} < L \x - y\ , 



with L :— \\3>(-,t-)\\ < oo, where boundedness follows from Lemma 2.2 and the 
compactness of I. □ 
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Note that|Lcmma 2.2| and |Lemma 2.3| did not use the Lipschitz continuity assumption 



in Definition 2.1 Even stronger, the continuity of x i— > (t i— > tp(t,t_)x) relies only on 
the continuity of tp as a function defined on I x I x R™ and mapping to R™. Although 
seemingly simple, |Lcmma 2.3 indicates an important difference in asymptotic and 



finite-time analysis: on unbounded (to the right) time-sets the continuous dependence 
of whole trajectories on the initial value corresponds to the definition of stability in 
the sense of Lyapunov, which is a nontrivial feature of certain trajectories. In the 
finite-time case, as |Lemma 2 .3| shows, it holds under very general assumptions. 

2.2 Logarithmic Difference Quotient 

The following concept will help us to introduce some of the forthcoming notions and 
to present the theory in an elegant and coherent way. 

Definition 2.4 (Logarithmic difference quotient). We define 

ln/(t)-ln/(.s) 



Ai: C(I,R >0 ) ->C(#ixi,E), / >-> (i, a) ^ 



t-s 



and we call Aj(/)(t, s) the logarithmic difference quotient of f at t and s for / G 
C(I, R>o) and t,s €^ixi- For notational convenience we will write A for Aj when 
there is no risk of confusion. 

Remark 2.5. Note that due to compactness of I any function / € C(I, R>o) is uni- 
formly continuous and /[I] C [a,b] with < a < b, i.e. / is bounded above and 
bounded away from zero. Furthermore, ln| ^ a ^ is Lipschitz continuous with Lipschitz 

constant - and bounded. Suppose / is additionally Lipschitz continuous, then 



|su P {A(/)(M); (M)e^ixi}|,|inf{A(/)(M); (*,«) e^ixi}| < 



I/Lip 



Furthermore, we have that A(-)(t, s) is continuous for any (t, s) S^ixi- Summarizing, 
A considered as 

A: C(I, R>o)x j^ixi -» R, (/, t, s) ^IRzMM 

t — s 

is continuous in the first argument and jointly continuous in the last two arguments, 
but, in general, not jointly continuous in all three arguments. However, for a linear 
process $ € CP(I, R n ) one can show that the family of functions (A(|$(-, t-)x\)) x ^s 
is equicontinuous. This can be calculated directly or, alternatively, can be deduced 
from lLemma 2.3l and the Theorem of Arzela-Ascoli. 

When applied to a linear process we recover joint continuity for the logarithmic 
difference quotient. 

Lemma 2.6. Let $ e CP (I, R"). Then 

^"\{0}x ^nxi^R, (x,t,s) i-> A(\<S>(-,t-)x\)(t,s), 

is continuous. 
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Note that the function defined in |Lemma 2.6| is well-defined, since we assume the 
linear process to be invertible (no trajectory starting in R™ \ {0} at t_ attains zero on 
I). 

Proof. Let x £ R™ \ {0}, (t, s) G^ixi and e £ R>o- By equicontinuity with respect to 
the initial value there exists Si £ R>o such that for all y £ R™ \ {0} 

max {|*-*'|,| S => \A(\$(-,t-)y\)(t >S )-A(\<Z>(-,t-)y\)(t',s')\<^. 



By Lemma 2.3 and hence continuity of R™ \ {0} 3 y H> A(|<I>(-, t_)y\)(t, s) there exists 
S 2 £ R> such that for all x' € R" \ {0} 

\x-x'\<5 2 \A(M.,t_)x\)(t,s)-A(\$(.,t-)x'\){t,s)\<~. 

Combining these two estimates, we obtain for (t' , s') G^ixi, x' £ R™ \ {0} 

\A(\^.,t-)x\)(t,s)-A(\$(.,t-)x>\)(t',s')\<£, 

whenever max{|t — t'\ , \s — s'\ , |a; — < min{<5i, £2}. □ 

Remark 2.7. The logarithmic difference quotient can be regarded as a finite-time ex- 
ponential growth rate as introduced in [9]. With the notation used there, we find 
that 

\^- s \sMs,t-)x) = Am-,t_)x\){t, S ). 

[Definition 2.4| and |Lemma 2.6| will have a major impact on the following theory 
and allow a development of the known finite-time spectral theory, starting from the 
notion of the logarithmic difference quotient alone. Another useful observation is the 



following, which holds obviously by linearity of s) and Definition 2.4 



Lemma 2.8. Let $ e CP (I, R"). Then for any iel"\ {0} and X £ R\ {0} one has 

A(M;t-)(\x)\) = A(M;t-)x\). 

As a consequence, we can equally well define A(|<I>(-, on the quotient set P n , 

the (real) projective space obtained by identifying vectors x, y £ R™ \ {0} if they are 
nontrivial multiples of each other. The projective space P n_1 can be identified to the 
set of one-dimensional subspaces, or lines through the origin, in R™. A generalization 
of that concept is the Grassmann manifold Gr(fc,R") of /c-dimensional subspaces of 
R", which will turn out to be very useful. Let us briefly recall the construction and 
its basic properties (following essentially jTTJ [T] ) . 

For any k £ {1, ...,n} one can introduce the Grassmann manifold Gr(fc,R n ) as 
follows: consider the set of all orthonormal k- frames in R™, which is also referred to as 
the Stiefel manifold St(fc,R"), i.e. fc-tuples of orthonormal vectors in R™ represented 
by n x fc-matrices A such that A T A = id. It is well-known that the Stiefel manifold 
is a compact subset of R nxfe . The Grassmann manifold can then be obtained as 
the image of the function ir: St(fc,R") — > Gr(fc,R n ) mapping A = (ai • • • ak) to 
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span {ai, . . . , a^}. Endowing Gr(k, R n ) with the final topology with respect to tt turns 
it into a compact set. As is proved in [T7], the so-called gap metric 8: Gr(fc,R") 2 — > 
R>o, Q(L,M) = \\ttl — 7rj\4-||, where ||-|j denotes the operator norm and irx denotes 
the orthogonal projection onto a subspace X, induces a topology that coincides with 
the final topology of tt. Furthermore, the same topology is induced by the Hausdorff 
distance between the intersections of, respectively, the subspaces L and M with the 
unit sphere (see, for instance, [THl Chapter 13]). As a consequence, X n "~ i ' 00 > X for 
X e Gr(k, R") and (A„) nS]N G (Gr(fc, R n ))® implies that for any x e X there exists 
a sequence (x n )n&i € with a;„ G X n for each nel such that x n "^°°> x with 
respect to any norm on R" by the equivalence of norms in finite-dimensional normcd 
vector spaces. 



As a consequence of Lemma 2.8 we find that x n- A(|$(-, t_)x\)(t, s) is actually 
uniformly continuous for any (t, s) G^ixi- The following simple observation associates 
the notions introduced in the present work to notions with the same nomenclature 
appearing in the references HH [13] . 

Lemma 2.9. Let f 6 C(I, R>o)- Then the following statements are equivalent: 

(i) There exists 8 £ R>o such that for any t,s G I, t > s, one /ias /(i) < 
e-*(*-«)/(a). 

(ttj There exists 8 € R>o such that t H >• e St f(t) is decreasing, 
(in) sup{A(/)(i,s); (t,«) G^ixi} < 0. 

Moreover, if I is an interval and f is differentiable, then each of the above statements 
is equivalent to 

(iv) there exists 8 € R>o sucA £/iai /' < (5/ holds. 
Analogue statements hold for / e C(I, R >0 ) if inf {A (/)(*, s); (t,s) G^ixi} > 0. 



2.3 Growth Rates 

Due to lLemma 2.81 we can define 

R" \ {0} x ^ Ixl -+ R, (x, i, s) ^ A(|*(-, t_)z|)(*, a), 

equivalently by 

P" _1 x ^nxi R, {[x]~,t,s) i ^ A(|*(-,t_)s|)(t,«), 

and, moreover, by 

Gr(l, 1") x ^ Ixl -> R, (X, t, s) = (span {a;} , t, s) h> A(|*(-, i_)x|)(t, a). 

On the other hand, we can identify an equivalence class P™" 1 9 [a;]^, = X E Gr(l, R") 
with its single normalized representative lying in some half-space, for instance R>o x 
R"" 1 . By [H p. 198, Eq. (2.12)] the gap metric on Gr(l,R") and the metric induced 
by |-| are equivalent. 
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Definition 2.10 (1-dimensional growth rate). Making use of the above identifications 
we define 



A 1 , A : Gr(l,R") x £P{I,R n ) -> R, 
(X, $) = Qx]„,*) H- A r (X, $) := inf {A(|$(-, t_)as|)(*, s); (t, s) G^ixi} , 
(Y, $) = (|y]„, $) i ^ AV, * ) - sup {A(|$(-,t_)j/|)(i, s); (t, s) e^ lxl } , 

and call A I (X, <I>) and A (X, <fr), respectively, the lower and upper growth rate of X 
under $. For convenience, we drop the index I in case the time set is clear from the 
context. 

One readily verifies that for any $ £ £P(I,R n ) and X £ Gr(l,R") one has 
|A(X, $)| , A(X, $) < L1-L2 < oo, where Li e R is the Lipschitz constant of the 



logarithm restricted to some interval as in Remark 2.5 and L 2 £ R is the (global) Lip 



schitz constant from [Definition 2. 1| The boundedness of growth rates for nontrivial 
subspaces does not hold in general when we weaken the assumed Lipschitz continuity 
of linear processes to absolute continuity, see [7] . 

Remark 2.11. One can consider the 1-dimensional growth rates as finite-time analogues 
to Bohl exponents, see [5] and also [T2J p. 118, pp. 146-148] for the definition and a 
historical review. 

Proposition 2.12. Let $ £ CP (I, R"). then A(-,$) and A(-,<I>) are continuous and 
bounded. 

Proof. We prove the assertion for A(-,<£) since the argumentation is completely anal- 
ogous for A(-, $). To this end, we prove upper and lower semi-continuity separately. 

The upper semi-continuity follows from the identification of 1-dimcnsional subspaces 
with their unique normalized element lying in the half-plane, the continuity of x i— > 



((t,s) i — y A(|$(-,t_)a;|)), see Remark 2.5 and the well-known fact that the infimum 
function over families of continuous functions is upper semi-continuous. 

We prove lower semi-continuity of A (•,<£) by contradiction, assuming that there 
exist some subspace X £ Gr(l,R") and a sequence (Xi) i& K £ Gr(l,R") w with 
lirrij^oo Q(Xi, X) = and a' := liminf^oo A(Xj, $) < A(X, $) =: a. Define rj := 
a ~ a £ R>o- Then there exists a subsequence (Xi,)j^ such that for each j £ IN we 
have ApQ., $) < a' + r) = a — n. Let (xj)j & n £ S N be the associated sequence of 
characterizing elements, then we can find sequences (tj, s,-)j £ ]n £ (^ixiO^such that 

MM;t->Mtj^j)<a-V (2-2) 

for any j £ IN. By assumption, we have that Xj J ^°°> x £ S n X converges to the 
characterizing element of X, and due to the relative compactness of =/=i x i we may 
assume without loss of generality that 
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Our aim is to find (u,v) G^ixi such that A(|<&(-, t_)x\)(u, v) < a by inequality (2.2) 
and some continuity argument. To do this, we need to distinguish two cases: (1) 
t* ^ s*, i.e. (#%«*) e^ixi, and (2) t* = s* . 



(1) In this case, we can choose (it, v) = (i*,s*), since by Lemma 2.6 and inequality 
(2.2) we find that the limit still satisfies the inequality, i.e. 

A(m;t-)x\)(t*,s*)<a-r), 

which contradicts the minimality of a — X(X, <&). 

(2) First note that in this case, the limit points are not in the domain of the A-function 
as defined in |Lcmma 2.6| Nevertheless, t* is necessarily an accumulation point in 
I. The idea is to perturb t* without altering the logarithmic difference quotient 
too much. There exists a function jo : R>o — > IN such that for any 6 € H>o such 
that |tj (<5) — 1*\ < 8 due to the fact that t* is an accumulation point. By the 
convergence of (tj)j and (sj)j there exists a function j\ : IN x K, >0 — > IN such that 
for any j e N>j 1 (i (5/2),5/4) n °l ds 

s s s 

\*j - tjol + l«j - s *\ < \*J -t*\ + \t*-th\ + \sj-s*\< - + - + - = 5. 



By Lemma 2.6 we have on the one hand that there exists a function Si : R>o — > lR>o 
such that for any e € R>o and y € S the inequality \tj — tj \ + \sj — s*\ < Si(e) 
implies 

\A(M-,t_)y\)( tj , S3 ) - A(\$(-,t_)y\)(t JO , S *)\ < e. 

On the other hand, by the same lemma, there exists a function j 2 : H>o — ► IN such 
that for all j £ 1N >J2 ( £ ) 

|A(|$(-,t_)^|)(t io , S *)-A(I^O.*-W)(*io,s*)|< e . 
Combining the above estimates, we obtain 

|A(|$(.,t_)a;,|)(t J , Sj )-A(|cI>(.,^) a ;|)fe , S *)| 

< |A(|*(.,t_) aJi |)(t ij ^)-A(|*(. J t_) a : J -|)(i i0 ,a*)l + 

+ |A(|*(.,i_)a: i |)(i i0J O-A(|$(- J t-)x|)(i i0J OI 

< 1 + 1 < 1 
~ 4 4 - 2' 

for all j > N := max{j ( S -^),j 1 (j (S 1 (ri/4),S 1 (r ] /4)),j 2 (T 1 /4)y For j e TN >N 
we estimate 

Ad^C-,^)^!)^,^) < |A(|<&(-^_)x|)(i, ,^) - A(|3>(-, £ _)^|)( £j , Sj -)| + 

+ A(|$(- ) t_)a: J |)(* J) a J -) 

?7 T] 
<a-T}+-=a--, 

which contradicts the minimality of a — X(X, <I>). 
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Hence, we have 

liminf \(Xi,®) > A(X,$) 

i— >oo 

and the assertion is proved. □ 

For later robustness investigations the following (semi-)metric and the induced topol- 
ogy on CV(I, R") will turn out to be very helpful. 

Definition 2.13 (Metric on CT(I,R n )). We define di,df. CV(I,R n ) 2 -> R> by 
4: ($,*)>->• sup (max {| A(X, $) - \(X,V)\ , | A(X, $) - A(X,*)|}|, 

X£Gr(l,E") ^ ' 



dj : ($, *) i ^ max <^ sup || (*(., i_) - ¥(., *_))»!!„ , di($, ¥) 

Obviously, dj is a semimetric, i.e. it satisfies non-negativity, symmetry and the triangle 
inequality (but not definiteness), and dj_ is a metric on £7 7 (I, R"). 

Note that, in general, di cannot be extended to a proper metric on the set of abso- 
lutely continuous linear processes due to possibly unbounded growth rates. Neverthe- 
less, the (semi-)metric can be used to define open balls around absolutely continuous 
linear processes. 

In the following, we endow CP (I., R") with the topology induced by dj. We extend 
the definition of 1-dimensional growth rates to higher-dimensional subspaces. Since 
both definitions coincide for 1-dimensional subspaces, we use the same symbol without 
risking confusion. With Lemma 2.9 one easily establishes equivalence to the growth 
rates as defined in [7j [14j [13] . 

Definition 2.14 (Growth rate, cf. 0Q1CE3]). We define 

n 

A 1 : (J Gr(fc,R n ) x CV(I,R n ) -> R, {X, $) i-> min A I ([Gr(l, X)}, 
fc=i 

n 

A 1 : (J Gr(fc,R") x £P(I,R n ) -> R, (X, $) i-> max A I ([Gr(l, X)], $), 
fc=i 

and call A E (X, and A E (X, $), respectively, the lower and upper growth rate of X 
under <f>. For convenience, we drop the index I in case the time set is clear from the 
context. 

Note that the above definition is well-defined by the continuity of A(-, $) and A(-, $) 
and the compactness of Gr(l, X). We extend the definition naturally by A({0} , $) = oo 
and A({0} , <£>) = — oo for any $ € CP(I, R n ). The following simple observations follow 
directly from the definition. 

Lemma 2.15 (cf. [14, Remark 6]). Let $ e £"P(I,R n ) and X, Y C R" be two non- 
trivial subspaces. Then 
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(i) X C Y implies X(Y, $) < X(X, $) and X(Y, $) > A(X, $); 

fuj inr / {0} implies X(X, $) < A(F, $) and X(Y, $) < A(X, $). 

Proof, (i) This is clear since, with respect to Y, infimum and supremum are taken over 
a larger set, respectively. 

(ii) Let x G (X n Y) \ {0}, then by definition we have for any (t, s) G^i x i 

A(X,$) < A(|$(-,t_)z|)(M) < A(y,$) 

and 

A(y,$) < A(|$(-,t_):c|)(M) < A(X,$). □ 

For the general growth rates restricted to some Gr(fc,H n ) we obtain continuity as 
well. This extends [TU Theorem 9] and will be used in [Proposition 4.6 to generalize 
QH Theorem 14]. 

Proposition 2.16. Let $ G CV(I, R n ). Then for each k G {1, . . . , n}, the restrictions 

A(-> $ )| G r(fc,]Rn) 0nd I (-' $ )lGr(fe,K") 

are continuous and bounded. 

Proof. Let I,ye Gr(/c,R") and e G H>o- We prove the assertion only for the lower 
growth rate, since the proof is analogous for the upper one. Without loss of generality 
we assume that X(X, $) < A(y,$). Let U := argmin A([Gr(l, X)], $) G Gr(l,X) and 
weSbc its characterizing element. Choose (t, s) G^ixi such that 

A(|$(-, t_)«|)(t, s) - A(C7, $) = A(|*(., t-)u\)(t, s) - A(X, <D) < e/2. (2.3) 

By the uniform continuity of x H> A(|$(-,t_)ar|)(t, s) there exists 5 G R>o such that 

|A(|$(. ) t_)a:|)(*,5)-A(|$(. ) *_) I ,|)(*, fl )|<e/2 (2.4) 

whenever \x - y\ < 8, x, y G 5. If 8(X, F) < 6/2 then by [Ml P- 198, Eq. (2.12)] there 
exists V G Gr(l, Y) such that for its characterizing element v E V P\ S the estimate 
|u — v\ < S holds. Combining Eqs. ( |2.3[ ) and (2.4), we obtain 

A(Y, <D) - ApT, $) < A(|*(-,i_)i;|)(t, «) - A(C/, $) 

<|A(|$(.,t_)«|)(t, fl )-A(|$(. ) t_)u|)(t, a )| + 
+ |A(|$(-,i_)u|)(t,s)-A([/,$)| <e, 

which proves the continuity. Boundedness follows from the compactness of Gr(fc, R"). 

□ 

By construction of the topology we obtain the continuous dependence of growth 
rates on the linear process. 

Corollary 2.17. Let X C R™ be a subspace. Then X(X, •), X(X, •) : £P(I,R n ) -> 
RU {—00,00} are Lipschitz continuous with Lipschitz constant 1. 
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2.4 Extremal Growth Rates 



The following notion is connected to subspaces of different dimensions with optimal 
growth rates and it will play an important role in the linear spectral theory as well as 
in robustness issues. 

Definition 2.18 (Extremal fc-growth rates, cf. HU [T3]). For fc G {1, ...,n} we 
define 

Xl: CP (I, R n ) -> R, $H^maxA I ([Gr(fc,R n )],$), 
X\: CP (I, R n ) -S-R, $ i-> minA I ([Gr(fc,R n )],$), 

and call A^(<i>) and A fc ( < &), respectively, the maximal lower and minimal upper k-growth 
rate of For convenience, we drop the index I in case the time set is clear from the 
context. We extend the above definition naturally by Aq($) = oo and Ao( < f ) ) = — oo. 

Remark 2.19. Note that the extremal fc-growth rate functions are well-defined due to 
the continuity from Proposition 2.16 and the compactness of (Gr(fc, R n ), O). That 
means in particular that for any <f> G CP (I, R") and k G {0, . . . , n} there exist X,Y E 
Gr(fc, R") such that X(X, $) = A fc (4>) and A(Y, $) = A fe ($). We refer to such subspaces 
as extremal subspaces. Of course, in general, extremal subspaces need not be unique. 

Lemma 2.20 (cf. 7, Remark 11], [14, Remark 9]). For any $ G CP(I,R n ), the 
extremal growth rates are ordered and nested as follows: 

- oo = A (<i>) < A„($) < A x ($) < A„_ 1 ($) < . . . 

• • . < A n _i($) < Ai($) < A„($) < Ao($) = oo. 

Proof. The ordering properties follow directly from Lemma 2.15[ i) and the nesting 
property from Lemma 2.15 ii). □ 

Another observation is the continuity of the extremal growth rate functions, which 
easily follows from |Corollary 2.17[ 

Lemma 2.21. For any k G {1, . . . , n} the extremal k-growth rate functions A fe , Afc are 
Lipschitz continuous with Lipschitz constant 1. 

Note that [13l Theorem 20] is essentially the e-5-notation of the continuity estab- 
lished in |Lcmma 2.21| for the special case that I is finite. 

More interesting is the following convergence result with respect to the time-set. 
We denote by g?h the Hausdorff metric on the space of compact subsets of M C R, 
denoted by K.{M). The proof in [TJ] applies also to our more general situation here. 



Lemma 2.22 ([13 Theorem 17]). Let $ G £P(I,R"). Then for each k G {!,..., n} 



and any sequence JT = (JJi) i(F]N G /C(I) W of compact subsets o/I with dn(I, J;) t ~ > °°> 



one has 



lim 

i— f C30 



A^^k^-AU*) =0, and lim A fc '($l 2 )-A,($) 



= 0. 
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3 Spectral Theory for Linear Finite-Time Processes 



This section is devoted to the development of a spectral theory for linear finite-time 
processes, which is in the spirit of [301 [3S1 [7J. We start with a notion based on a 
certain dynamical dichotomy. We denote by P(R") the set of (linear) projections on 
II" . Throughout this section let $ G CP(f, R") denote a linear process on I. 



3.1 Exponential Monotonicity Dichotomy 

Definition 3.1 (Exponential monotonicity dichotomy). $ admits an exponential mo- 
notonicity dichotomy (EMD) on I (w.r.t. \-\) if there exists k G {0, . . . , n} such that 

A fe ($)<0<A„_ fe (<i>). (3.1) 

For brevity, we sometimes call a linear process $ on I admitting an EMD also (finite- 
time) hyperbolic. We call a hyperbolic linear process $ on I (finite-time) attractive if 
k = n and (finite-time) repulsive if k = 0. 

In the next lemma, we characterize an EMD on I. 

Lemma 3.2 (cf. [Jj). The following statements are equivalent: 

(i) $ admits an EMD on I with k G {0, . . . , n}. 

(ii) There exist subspaces X G Gr(fc, R n ) and Y G Gr(n-k, R n ) such that X(X, $) < 

o<a(y,$). 

(iii) There exists a projection Q G P(R") with rkQ = k such that A(imQ,$) < < 
A(kerQ,$). 

(iv) There exists a projection Q G P(R") and constants a, f3 G R>o such that for any 
t, s G I, t > s, and x G imQ, y G kerQ one /las 

|$(M-)z| < e~^ s ) |<&(M-H , 
|1>(M-)!/|>e WM |$( S ,(-)i 



Proof, (i) (ii): By Remark 2.19 there exist subspaces X, Y of appropriate dimension 
satisfying A(X, $) = Xggj < < A fe ($) = A(F, $). 
(ii) (i): By Definition 2.18 we obtain 

A fc (<&) < X(X, $) < < A(F, $) < A„_ fe (<i>). 

(ii) <S=> (iii) : The subspaces X and Y define a projection by im Q := X and ker Q := Y 
and vice versa. 

□ 



(iii) <S4> (iv): This follows directly from Lemma 2.9 
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Clearly, the subspaces / projections mentioned in Lemma 3.2 give rise to so-called 
invariant projectors, for instance P(-) = $(•, i_)Q$(t_, •), with corresponding proper- 
ties, see [7J H HI M EH- 



As can be seen easier from Lemma 3.2 Definition 3.1 includes some other (hyperbol- 
icity) notions as special cases. In the finite-time context so far only solution operators 
$ G CP (I, R") of linear differential equations on a compact time- interval I = [r, r+T], 
t G R, T e R> , i.e. 

x = A(i)a;, 

where usually A S C(I, L(R™)), have been considered. In this case, $ regarded as a 
linear process satisfies even continuous differentiability where we required just Lipschitz 
continuity in Definition 2.1 From Lemma 2.9 and Lemma 3. 2| follows directly that, 
setting I = [t,t + T], our definition of EMD coincides with the definition of M- 
hyperbolicity as defined in [6l [14] , (finite-time) hyperbolicity as defined in [TJ [5] and 
uniform hyperbolicity as in [20] . By setting I = J it is obvious that our definition 
of an EMD corresponds essentially to the finite-time hyperbolicity proposed in [15] , 
which generalizes in particular the nonhyperbolic (r, T)-dichotomy as suggested in 28J. 



From Lemma 2.9| and taking into account the monotonicity preserving property of the 
logarithm one easily concludes that nonhyperbolic (r, T)-dichotomy is equivalent to 
Definition 3.1 with I = {r, r + T}. Another finite-time hyperbolicity notion of the 



EMD-type which is based only on the start and end time-point can be found in the 
finite-time Lyapunov exponent approach, see for instance [211 13 1] . As is shown in [13] 
this approach is closely related to our finite-time hyperbolicity notion. 



k G 



By the nesting property of the extremal growth rates, see Lemma 2.15 we have that 
{0, . . . , 7i } in Definition 3.1 is uniquely defined. Of course, this does not imply 

are unique. Consequently, 



that the subspaces / projections mentioned in Lemma 3.2 



the notion of EMD is well-defined, at least up to rank of the projection, in contrast to 
the definition of a finite-time exponential dichotomy. 

In the following, we investigate for given $ G CP (I, R") the associated family of 
linear processes (<& 7 ) 7e ]a G CP(1, R n ) E which is defined for any 7 G R and s,t € I by: 



$ 7 (t,s) 



:= c-^- s 



(3.3) 



The motivation to study $ 7 comes from the fact that they arise naturally in the study 
of linear differential equations 

x = A{t)x, 

where, for instance, A G C(I, L(R")), I an interval, and the corresponding shifted 
linear differential equations for 7 G R 

i = (A(t)-7id B »)x. (3.4) 

Clearly, we have by definition that for any x G R™ and (t, s) G^ixi holds 

As a consequence of this observation, we obtain the following results. Let $ G 
CP(l,n n ) and 7 G R. $ 7 admits an EMD on I with projection Q G P(R") if and 
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only if A(imQ, $) < 7 < A(kcrQ, <£>). Furthermore, the extremal growth rates behave 
as the logarithmic difference quotient under the exponential weight, i.e. 

A fe ($ 7 )+7 = A fc (*) and A fc ($ 7 )+7 = A fc ($). 

Consequently, the following characterization of EMD for weighted linear processes 
holds. 

Corollary 3.3. $ 7 admits an EMD on I with fee {0, . . . , n} if and only if A^($) < 

7 < A„_ fc ($). 

Since the semimetric dj depends only on the growth rates, we get the following 
normalization property: di(<I>, <£> 7 ) = 7. 

3.2 Spectrum of Linear Finite-Time Processes 

Next, we introduce a finite-time spectral notion which is based on the EMD. 

Definition 3.4 (Finite-time dichotomy spectrum, cf. [251 HI US] ) • We define 

S 1 : CP (I, R n ) -> 2 B , $ i-> {7 G R; $ 7 does not admit an EMD on 1} , 

and call E I (<i>) the (finite-time dichotomy) spectrum and := R\£ I ($) the (finite- 
time) resolvent set 0/$, respectively. 



Remark 3.5. By |Definition 3.4 it is clear that for two compact sets J, I C R with JCI 
and $ G CV(I,R n ) we have that £ J ($|j, 2 ) C E I (<i>). 

Theorem 3.6 (Spectral Theorem, [Ml Theorem 10], cf. also [7} Hi]). Denote 

{i ,...,i d } := {j G {0,. .. ,n}; Aj($) < A„_,($)} , i fc < i fc+ i, fc G {0, . . . , d - 1} . 
Then the spectrum 0/ $ is the union of d disjoint compact intervals, i.e. 

and we call [A„_ ifc i ($), Ai fc ($)] t/ie A:-th spectral interval. 
Proof. This follows directly from |Corollary 3.3| by the fact that 

n d 

p\$>) = U(^( $ ).A„_ fc ($)) = \j(Km^ n -iM)> 

k=0 k=0 

since (A&(<I>), A n _ fe (<&)) = for G {0, . . . ,n} \ {io, ■ ■ ■ ,id}, an d the nesting property 
of the extremal growth rates stated in |Lcmma 2.20| □ 
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The assertion of |Thcorcm 3.6| remains basically true if we require only absolute 
continuity of $. In this case, however, the left-most and right-most spectral intervals 



may be unbo unded, see al so the remark after Definition 2.10 and [7J Theorem 17]. 



Theorem 3.6 together with Lemma 2.20 implies that is non-empty for any $ € 

cp{i,n n ). 

Due to the simple interval-structure of the spectrum where the endpoints of the 
intervals are the extremal growth rates, we obtain the new result of continuous depen- 



dence of E 1 on the linear process from Lemma 2.21 



Proposition 3.7. The spectrum function S : (£"P(I, M n ), dj) — > (/C(Il), g?h) is Lips- 
chitz continuous with Lipschitz constant 1. 

Additionally, we obtain a kind of continuity result of the spectrum with respect to 
the time set, corresponding to |Lcmma 2.22| 

Corollary 3.8 (Q31 Theorem 17]). For any (Ji) ielN G K,{tf one has 
lim d H (I,JFt) =0 => lira d H (E I ($), E J *($| 2 )) = 0. 
Proof. This follows from |Theorem 3.6| and |Lemma 2.22] □ 

3.3 Robustness of Hyperbolicity 

A desired property of hyperbolicity is its robustness under perturbations, which we 
establish in a very natural way and for the first time in such a general setting. 

Theorem 3.9 (Robustness of EMD). Let $ admit an EMD on I with fee {0, . . . , n}. 
Then any * G CV(I,R n ) with < min {-A fe ($), A„_ fc ($)} admits an EMD 

on I with k (and the same extremal subspaces). 

Proof. This is a simple consequence of |Lemma 2.21 and Definition 2.13 i.e. the fact 



that the extremal growth rates of $ and ^ differ at most as much as $ and ^ do. 
By |Corollary 2T7l the same holds for the subspaces realizing the extremal growth 
rates. □ 

Corollary 3.10. For given I the set 

{$ E CV(I,U n ); $ admits an EMD on 1} 

is open in £"P(1, 11" ) with respect to the topology induced by dj. 

Next, we are going to show that the estimate in |Theorem 3.9| is sharp and hence 
gives the maximal perturbation bound for EMD persistence. This can be interpreted 
as a hyperbolicity radius or, in case <& is an attractive linear process, as a stability 
radius. 

Theorem 3.11 (Hyperbolicity radius). Let $ admit an EMD on I. Then 

g:= dist({0},E I ($)) e H>o 

is the largest number such that any ^ € £P(I, R") with rfi($,\I>) < g admits an EMD 
on I. 



1G 



Proof. The fact that g is a number with the asserted property is clear by |Theorem 3.9| 
thus it remains to show that it is the largest one. This is indeed clear by the normal- 
ization property of dj, i.e. dj($, $ e ) = g and that <£> e does not admit an EMD on I. 
The latter is due to the fact that g e 5? (<£>). □ 



4 Linearization of Finite-Time Processes 

In this section we study C^-processes ip on I and their linearization $ along fixed 
trajectories. In particular, we are interested in local implications of finite-time hyper- 
bolicity of $ on the the original process ip. Throughout this section let ip € V(I, R") 
denote a C 1 -process on I. 



4.1 Linearization and Hyperbolicity of Finite-Time Processes 

Motivated by the classical theory, we introduce the following notion. 
Definition 4.1 (Linearization). We define for any x € R™ the following function 

$(t_,x): I x I -> L(R"), (t,s) i y $( f _^)(t,s) := d 2 tp(t, s, tp(s, t_, x)). 
We call *&(t_ jX ) the linearization of ip along (p(-,t—,x). 

Lemma 4.2. For any x G R" the function <&(t_ jX ) as defined in Definition ^.1 is a 
linear process on I. 

Proof. The cocycle properties including the invcrtibility are easily checked with the 
cocycle properties of <p and the chain rule of differentiation. Lipschitz continuity of 
&(t-,x) holds by definition. □ 

Definition 4.3 (Finite-time hyperbolicity, attraction and repulsion). Let xq € R n and 
®(t-,x ) ^ CV(I, R") be the linearization of >p along ip(-,t„,x ). We call p(-,t^ 7 x ) 
(finite-time) hyperbolic if *&(t_ jX o) admits an EMD on I. We call (p(-,t—,Xo) (finite- 
time) attractive / repulsive if < &(t_, a;o ) is attractive / repulsive. 

For an extensive study of finite-time attractivity with respect to the two-point time- 
set I = we refer to [IB]. To get a better geometrical understanding of the 
finite-time behavior of trajectories close to some reference trajectory we introduce the 
following cones, i.e. sets that are invariant under scalar multiplication. 

Definition 4.4 (Stable / unstable cone, cf. [14., p. 11]). We define 

V B : CV{I,R n ) -> 2 R '\ $H.(J{X e Gr(l,R n ); A(X, $) < 0} , 

V u : CPpL, R") -> 2 K " , $ h> |J {X € Gr(l, R n ); A(Z, $) > 0} , 

where we call I^ S (<I>) and V u ($>) the stable and unstable cone o/<I>, respectively. 



17 



Clearly, we can identify V s ($>) and F u ($) to subsets of Gr(l,R"). Formally, for a 
nontrivial subspace V C R™ we identify 

V = {J{ X( ^ Gr(l,R n ); X C V} ^ {X e Gi(l,R n ); X C V} = Gr(l, V), (4.1) 

and in that sense we can speak of a subspace being compact. By the continuity of 



the growth rate functions for fixed $ (see Proposition 2.16), we obtain the following 
result. 

Lemma 4.5. For any <I> G £P(I,R") ifte stable and unstable cone of $ are open 
sufesete o/Gr(l,R"). 

The next proposition is a generalization of [141 Theorem 14], showing that a differ- 
entiability assumption is in fact unnecessary. 

Proposition 4.6 (cf. [H Theorem 14]). Let $ G CV(I,R n ). 

(i) //$ admits an EMD on I with Q G P(R"), then 

imQ C F B ($), and ker Q C F u ($). 

(ii) Suppose there exist subspaces U\ C V s ($) and J7 2 C V™($) wi/j/i Z7i © f7 2 = R"- 
T"/ien $ admits an EMD with projection Q defined by imQ — U\ and ker Q = U2 ■ 



Proof, (i): This follows directly from Definition 3.1 and Definition 4.4 



(ii): We need to show that \(Ui,<&) < < \{U 2 , $). This is clear by the continuity of 
the 1-dimensional growth rate functions from [Proposition 2.12 and the compactness 



of Ui and U 2 . □ 

Remark 4.7. From the proof of part (ii) we can read off directly that the supremum 
(infimum) over upper (lower) growth rates with respect to an arbitrary compact sub- 
set of F s (<&) (y u ($)), considered as subsets of Gr(l,R n ), gives a negative (positive) 
number. 

In [24 for two-dimensional Hamiltonian ODEs, in [20 for general two-dimensional 
ODEs, in [21 for general three-dimensional ODEs and in [7 for n-dimensional ODEs 
stable and unstable manifolds are introduced as manifolds in the extended phase space, 
depending on some non-unique extension of a given differential equation on some com- 
pact time-interval to the whole real line. In these works, for any extension the stable 
and unstable manifolds have indeed a C 1 -manifold structure. In [TSl Definition 35] for 
two-dimensional ODEs so-called stable and unstable manifolds are introduced, which 
do not have a manifold structure. However, these objects are defined in an "intrin- 
sic" way, by requiring some decay and growth property of a solution on the compact 
time- interval with respect to a reference solution. In 18, Definition 3.1] domains of 
attraction are introduced intrinsically with a decay requirement with respect to the 
two-point time-set I = {/j_,t + }. We take this as a motivation for the next definition. 
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Definition 4.8 (Domains of finite-time attraction / repulsion). Let x e R™. Then 
we define 

W* :=\xeR n \{x }; sup {A(\<p{-, i_, x) - <p{; t_, a; )|)(t, *)} < \ U {z } , 
[ (t,«)e^ixi J 

^ u := | a; e R™ \ {x } ; inf {A(|^(-, t_, x) - <p(., i_, x )\)(t, s)} > oj U {x } , 

and call W| and W^ the domains of (finite-time) attraction and repulsion with 
respect to i_, #0), respectively. 

Remark 4.9. We choose to call VK! and domains of attraction and repulsion, 
respectively, to emphasize the bare set structure and to avoid terms like manifold or 
cone. 

It is easy to see that under a time-dependent coordinate shift (t, x) n- (t, x — 
ip(t,t-,Xo)) the linearization 3>(t_ lXo ) (along ip(-,t-,Xo) in the original coordinates) 
coincides with the linearization $( t _.o) (along E x {0} in the transformed coordinates). 
Without loss of generality we assume the reference solution to be the zero solution for 
the rest of this section. 

As a next step, we prove that hyperbolic trajectories have, under some condition on 
the approximation quality of the process by the linearization, non-empty domains of 
attraction and repulsion. We are going to show that locally cones and domains look 
very similarly, which is a result of local persistence and we adapt the reasoning that 
led to the robustness of EMD to the nonlinear case. To this end, we first introduce 
analogues to the 1-dimensional growth rates as follows: 

M: R"\{0} xP(I,1") ->R, {x,<p)i-> inf A(\<p(-,t-,x)\)(t, s), (4.2) 

(M)6#ixl 

ju: R™\{0} x V(I,R n ) -> R, {x,<p)^ sup A(\ip(-,t-,x)\)(t,s). (4.3) 

(M)e#ixi 

Note that for $ E CV(I, R") and iel"\ {0} we have fi(x, $) = A(span {x} , $) and 
~p(x, 4>) = A(span {x} , <£>). Based on /i and JL we introduce a measure of approximation 
of the C 1 -process ip by the linearization $ along the zero reference trajectory 

m: R> -> R>o, (4.4a) 

'0, 77 = 0, 

v Sup x£B[0;7)] \ {0} max{|/i(a;,^) -/x(a;,$)| , \ji(x, ip) - ~p(x, $)|} , otherwise. 

(4.4b) 



With this notation at hand, the domains of attraction and repulsion of the zero refer- 
ence trajectory take the simple form 

:= {x e R" \ {0} ; Ji(x, <p) < 0} U {0} , 
I^ U := {i£E"\ {0} ; fi(x, ip) > 0} U {0} , 
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from which the similarity to the stable and unstable cones of [Definition 4. 4| becomes 
already visible. 

Next we give a sufficient condition for m to be continuous in 0. 

Lemma 4.10. Let $ £ CV(1, H n ) be the linearization of ip along ip(-,t-,0) =0. If I 
is finite then the function m as in Eq. (4.4) is continuous in 0. 



Proof. By the inverse triangle inequality it suffices to prove 

D(t,s,x) :=\A(\<p(;t-,x)\)(t,s) - A(m-,t-)x\)(t,s)\ 

uniformly in [t, s) S^ixi- Since I is finite we have M '■— 
Observe that 



> 0, 



M I* - *l > o. 



R(t,x) := (p(t,t-,x) - d 2 cp(t,t-,0)x = <p(t,t-,x) - <S>(t,t-)x 
is continuous and 

\R(t,x)\ > 0, uniformly in tel. 

We estimate 



(4.5) 



D{t, s, x) 



\a.\ip{t,t-,x)\ - \n.\<p(s,t-,x)\ In t-)x\ - In |$(s, t-)x\ 



< 



1 



M 
1 

~ M 



In 



t - s 
\<f{t,t-,x)\ 



t- s 



\9(t,t-)x\ 



111 



\cp(s,t_)x\ 



\*(8,t-,x)\ 



hi 



Mt,t-)x\ 



111 



|$(s,t_)a; + ii(s,x)| 



|$(s,t_)a:| 



>-(lnl + lnf) 



0, 



uniformly in (t, s) S^ixi by Lemma 2.2 and Eq. (4.5). 



□ 



Note that in the previous lemma we did not impose any extra regularity conditions 
neither on the norm nor on the process. The next lemma gives sufficient conditions 
for the ODE case, which requires both additional regularity of the norm and of the 
process. We state and prove it for the Euclidean norm, making use of the following 
facts: the Euclidean norm is continuously differentiable on W 1 \ {0}, and the modulus 
of continuity of its derivative when restricted to a compact domain D not containing 
the origin is u(t) — t/a, where a := min|[D]|. The following result can be clearly 
transferred to other norms by requiring continuous differentiability and a certain be- 
havior of the modulus of continuity of the derivative close to the origin, which will 
become clear in the course of the proof. 

Lemma 4.11. Let I be a compact interval, \-\ denote the Euclidean norm, ip £ 
V(l, El n ) be a C 2 '-process on I and $ e £7 7 (I, R n ) be the linearization of <p along 
(p(-,t—,0) = 0. Suppose that 



(i) (p and $ are continuously differentiable in the first argument, and 
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(ii) for R(t, x) := cp(t, t_,x) — <£>(£, tJ)x, t G I, x G K™ ; one /ias £/ia£ i?(i, •) together 
with doR(t, •) is of class 0(|a;| 2 ) for \x\ — > uniformly in t G I. 



TTien £/ie function m as in Eq. (4.4) is continuous in 



Proof. As in the previous lemma, we show that D(t, s, x) — > as — > uniformly 
in (t, s) G^ixi, this time by applying the mean value theorem to the continuously 
differentiable function 1 1— > In plfe^CT ■ Besides elementary calculations and estimates, 
the crucial ingredient of the proof is to show that 

\\\'\'m,t-)x + R(t,x))-\^m,t-)x)\\^^0, (4.6) 

uniformly in t G I. To show this, we first observe that there exist constants s, 5 G 
R>o such that \R(t,x)\ < e\x\ 2 < e5\x\ whenever |x| < S, due to the conver- 
gence assumption on R. Since <3? is invertible and I is compact, we have that a '■= 
min ||$( [I], £_)[£> ]| ,0 ■— max ||$([I], t_)|| > 0, where a is the absolute value closest to 
zero that a trajectory starting on the unit circle attains on I and, analogously, /3 is 
the largest such value. We may assume w.l.o.g. that S < a/e and hence a — eS > 0. 
Now choose n < 5, then for any |x| = rj we have that |$(t, f_)x| G [ar),/3rj\ and 
\tp(t,t-,x)\ G [(a — e8)ri,(3rj + en 2 ] =: I v . When restricted to the compact annulus 
B[0,f3n + en 2 ] \ B(0, (a ~ eS)n), the derivative of the Euclidean norm is uniformly 
continuous with a modulus of continuity of uj(t) = t/(a — s5)rf). On the other hand, 
on this annulus and all annuli constructed for smaller n we have the quadratic estimate 



on R, yielding uj(en ) = erf /{a — eS)n) and in turn proving Eq. (4.6). □ 



4.2 The Local Stable / Unstable Cone Theorem 

In the next two sections we introduce a novel, "intrinsic" approach to local stable 
and unstable cones and manifolds, which uses information about tp on I only and 
hence does not rely on classical asymptotic methods. The crucial assumption is the 
continuity of m in 0, which we have shown so far only for the case that I is finite. 
We define the two functions 

n: Gr(l,R n ) -> R> , X i-> inf {r G R>o; xeXnS, rx i W n s } , 

r): Gr(l,R") -> R> , X i-> inf {r G R>o! xeXnS, rx $ W£} . 

Theorem 4.12 (Local Stable / Unstable Cone Theorem). Let $ G CV(I,M n ) be the 



linearization of tp along cp(-,t-,0) = 0. If the function m as defined in Eq. (4.4) is 
continuous in then for any X, Y G Gr(l,R") with X C l/ s ($) and Y C V u ($) 
one has r)(X),f)(Y) > 0. Moreover, n and fj are bounded away from zero on compact 
subsets o/V rs ($) and V n (<$>), respectively. 



Proof. By Definition 4.4 we have for any X G Gr(l, R") with X C y s ($) that 
A(X, <£) < 0. By the continuity assumption on m there exists 8 G R>o such that 
771(77) < —X(X, <f>) for any 77 G (0, 6]. Then for any x £ B[Q, 5]D X we have Jl(x, tp) < 
and analogously the assertion for Y G Gr(l,R n ) with Y C U u ($). The second part 
follows from |Rcmark 4.7| and the same argument as applied before to single directions 
X G Gr(l,R n ). □ 



21 



By the same continuity argument as in |Theorem 4.12| we can find positive radii 
5 £ R> for the directions in the interior of V := It" \ (y s ($) U U u ($)) such that 
V n B[0, S] n W|,Fn B[0, 5] n W£ = 0. Roughly speaking, we find that stable and 
unstable cones of the linearization <£> and domains of attraction and repulsion of the 
process ip, respectively, are locally indistinguishable. Note that this is a pure continuity 
result and not an implication of hyperbolicity. For the ODE case in R 2 and stronger 
regularity assumptions a similar approximation result has been proved in [151 Theorem 
44]. 

As a special case we obtain the following result. 
Theorem 4.13 (Local Stable / Unstable Manifold Theorem). Suppose the assump- 



tions of \Theorem ^.12 are satisfied and let Q G P(R") be a projection such that 
ivaQ C V s ($) and kerQ C V n (<&). Then there exist neighborhoods U and V of the 
origin such that im Q n U C Wq and ker Q n V C Wq . Furthermore, for any t £ I i/ie 
following tangencies 

T ip(t,t-,[imQ]) =*(t,i_)[imO], and T ol p(t,t-, [kerQ]) = $(t,i_)[ker Q], (4.7) 

are satisfied. Consequently, </?([!], t_, [imQ n [/]) and y([I], t-, [kerQ n V]) can oe 
considered as finite-time local stable and unstable manifolds, respectively. 

Proof. Since imQ and kerQ are compact subsets of U S ( ( E>) and V ru ( < i>), respectively, 



Theorem 4.12 applies and the first part is proved. Furthermore, the tangencies (4.7) 



are easily verified with the definition of the linearization. □ 

Remark 4.14. We want to comment on some issues concerning |Thcorcm 4.13] 

1. |Theorem 4.13| holds for general compact I. So far, finite-time Local Stable Man- 
ifold Theorems have been proved only in the ODE case with I being a compact 
interval (see [2H[HJ[S], see a l so Remark 5~9^ 2). 



2. In the finite-time context, one can consider |Theorem 4.12| and |Thcorem 4.13| 
as robustness results as well. As we proved, locally, the stable and unstable 
cones (together with the subspaces that they contain) persist under nonlinear 
perturbations with vanishing first order terms. 

3. Despite the lack of structure for the domains of attraction and repulsion them- 
selves, we see that the maximal dimension of manifolds contained in these do- 
mains going through the origin corresponds to the indices of the EMD growth 
rates, i.e. to rank and deficiency of the EMD-projection, respectively. For short, 
domains of attraction and repulsion of hyperbolic solutions behave as desired in 
the sense of dimensionality. 

4. Note that by the assumption that ip be a C fc -process on I, k £ {1, 2}, we obtain 
directly that the extension of im Q n U and ker Q n U by <p to the extended state 
space I x R™ gives a C fe -manifold in each time- fiber {t} x R ra , tel. Evidently, 
the chart is given by <p(t,t-, ■). 
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5. The function m can be considered as a local measure of nonlinearity of ip, in 
the sense that m = if cp itself is linear and that m takes small values in case 
tp is only a small perturbation of a linear process. The more linear tp becomes, 
the larger we can choose the radius S £ R>o such that B(0,6) nimQ C Wq 
and B(0,5) n kerQ C Wq. In the "linear limit" we recover that imQ C Wq 
and kerQ C Wq. In this sense, we believe that our version of a local finite-time 
stable manifold theorem is the most natural one. On the other hand, for fixed 
nonlinearity, clearly m is an increasing function, i.e. the further away we go from 
the hyperbolic reference trajectory, the weaker we expect the exponential decay 
/ growth to be, until some point where the EMD-subspaces / cones leave the 
domain of attraction and repulsion, respectively. 

As an easy consequence of |Theorem 4.13| we obtain the following finite-time analogue 
of the classical theorem of linearized asymptotic stability. It is a generalization of [UJ 
Theorem 5.1] to (at least) arbitrary finite time-sets. 

Theorem 4.15 (Linearized Finite-Time Attraction / Repulsion). Let tp £ V(I, R") 
be a ^-process on I, x £ R", ip(-,t—,x) an attractive (repulsive) trajectory and m be 
continuous in 0. Then there exists a neighborhood U of x with U C WJ (U C W£, 
respectively) . 



4.3 Extensions of Cones and Domains and their Relationship 

Next we investigate the relationship between the cones V s ($) and V u (&) extended by 
the linearization $ to the extended state space on the one hand, i.e. 

V| := $(•, i_)[F s ($)] = {(t, $(t, t-)x) Six R"; (i, x) £ I x V S (<S>)} , 
V£ := $(•, t_)[V ll {<S>)} = {(t, $(t, t-)x) GlxR"; {t, x) £ I x V u {$)} , 

and the domains Wq and Wq extended by <p on the other hand, i.e. 

W S := v>(-,t_, [W S ]) = {(t,<p(t,t-,x)) £ Ix 1"; (t,x) £ I x W S } , 
W U := V (.,t_, [W U ]) = {(t, p(f, L.xjJelxE"; (M)elx^l- 

As usual, we denote by V|(t), V£(i), Wg(t) and W u (i), f £ I, the t-fibcr of the 
respective subsets of the extended state space. The next proposition states that in each 
i-fiber the extended stable and unstable cones are locally contained in the domains of 
attraction and repulsion, respectively. 



Theorem 4.16 (Relationship between Extensions). Suppose the assumptions of The- 
\orem ^.Tl| are satisfied. Define 

r7:IxGr(l,R")^R> , [t, Y) ^ inf {r £ R >0 ; y£YC\S, ry^W s (t)}, 
r):IxGr(l,R")^R> , (t, Y) n- inf {r £ R> ; y£YnS, ry<£W£(t)}. 

Then for each t £ I one has rj(t, -)\ vs ^yf](t, Olyift) > 
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Proof. We prove only rj(t, 



IV|(t) 



> 0, since the second assertion can be proved com- 



pletely analogously. Let t £ I, y £ V|(t) Pi 5. Since the proof is rather technical, for 
convenience, we sketch its idea first: By the invariance of V| under <!> it is clear that 
$(t-,t)y £ V|(i_) = y s ($). Now, consider a; r := <p(t-,t,ry), r £ (0,1]. To prove 
that ry £ Wo(i) for sufficiently small r, it is sufficient to show that a; r <G Wq = Wg(i-)- 
Since y s ($) is open by Lemma 4.5 our aim is to show that for sufficiently small r 
we have that x r is contained in a neighborhood of 1-dimcnsional subspaces around 
span {$(£_, t)y) which is a subset of V s (§). By Theorem 4.12 we then conclude that 
for sufficiently small r the vector x r is contained in the domain of attraction Wq. This 



proves the strict positivity of rj(t, •) 



V|(t) 



as claimed. 



Thus, it remains to show that x r £ Wg(i-) for sufficiently small r. By Lemma 4.5 
there exists 8 £ R>o such that 

B@(t-,t)y,9)c VUt-) = V s ($). 

Clearly, due to the positive homogeneity of the Banach space norm |-| on R™, the 
invariance of VJ(£_) under scalar multiplication and linearity of $(i_,t), we find that 
for all r £ (0,1] holds B($(t-,t)(ry),r9) C V|(t_). By expanding tp(t-,t,z) in 
z € R" \ {0} around we obtain 



\<p(t-,t,z)-$(t-,t)z\€o(\z\) for \z\ 



0. 



This is equivalent to the fact that for any e £ R>o there exists <S £ R>o such that for 
any z £ R" \ {0} with \z\ < S we have 

| V (t_,t, «)-*(*_, t)«| <e|z|. (4.8) 

In other words, for any e £ R>o, sufficiently small 5 £ R>o and \z\ < (5 we have 

tp(t_.t,z) £ B(®{t_,t)z,e\z\) C B{z,eS). 

and conse- 



In particular, choosing e = 0/2 we find that for 5 £ R>o from Eq. (4. 
quently for any ry with r €E [0,min{(5, 1}] holds 

x r = <p(t.,t,ry) £ B[<S>(t_,t)(ry),r9/2] C B($(t_,t)(n/),rfl) C V|(t_). 

Recall that B [span {$(t_, i)(ry)} , r6*/2] C Gr(fc,R") is closed and hence compact. 
Hence, by |Theorem 4.12| we know that for sufficiently small r £ R>o we have x r £ 



5 Applications 

In the following, we want to apply the above notions and results to ordinary differential 
equations, which we define without loss of generality globally for the sake of simplicity, 
i.e. 

x = f(t,x), (5.1) 



24 



where / G C 0,1 (I x R n ,R n ), J CI and (f(t, -))tei is uniformly Lipschitz continuous 
with Lipschitz constant Lf. Thus, Eq. ( |5.1| ) is well-posed and the solution operator 
tp is well-defined. It is well-known that ip satisfies the conditions of a C 1 -process. We 
consider the case that |-| is continuously differentiable on R™ \ {0}. In particular, this 
covers all norms induced by the Euclidean inner product and a symmetric positive 
definite matrix T G R™ xn as considered in (4[ [5] . We fix a solution ip(-,t-, x) : I — > D, 
x G D, and perform a time-dependent coordinate shift of the form (t 7 x) H t (t,x — 
(p(t,t-,x)) —. (t,y). Then in the new coordinates Eq. (5.1) takes the form 

y = d 1 f(t,<p(t,t-,x))y + g(t,y)=A(t)y + g(t,y), (5.2) 

where A := d x f{-, </?(•, t-, x)) G C(I,L(R")), g G C(I x R™,R") and 

g(t, v) = f(t,v + ip(t,t-, x)) - f(t,ip(t,t-,x)) - dif(t,<p(t,t-,x))v, tel, veU 

for t G I is the nonlinear term. By definition of the derivative we have 

g(t,x)/\x\^^0 

for any tel. In other words, for any t G I and e G R>o there exists S G (0, 1] such 
that the estimate sup {\g(t, x)\ / \x\ ; x G B(0, S) \ {0}} < e holds. Due to the uniform 
continuity of <?| IxB [ x y w e even obtain that for any s G R>o there exists S G (0,1] 
such that 

sup |^|^; xeB(0,S)\{0} : telj <e. (5.3) 

As usual, we call 

y = d 1 f(t,<p(t,t.,x))y (5.4) 



the linearization of (5.1 1 along (^(-, t_, x). It is well-known that the associated solution 
operator $ of (5.4), interpreted as a linear process, is the linearization of <p along 
(p(-,t—,x) and that $ is continuously differentiable in the first argument. Under the 
differentiability assumption on the norm the growth rates take the form (as introduced 
mH) 



A(X,$) =min( (l ^'/- )g ' ) ', ( * ) ; tel, x G X D S 
I \$(t,t-)x\ 

\(X,$) =max( %^M ; tel, i£ln5 



(5.5) 



for X G Gr(/c, R"), which can be seen by |Lcmma 2.9 From this representation and the 
chain rule (|$(-, t_)x\)'(t) = |-|' ($(t, t_)x)9 ($(t, t_)x) for f e I we can see that 
(dj) can be interpreted as some kind of C 1 (semi-)metric for linear solution operators 
on I. 

By classical techniques and Gronwall's lemma one easily establishes that linear right 
hand sides A G C(I, L(R™)) map continuously (with respect to c?i) to their unique 



solution operator $ G £"P(I, R"). From this observation together with Theorem 3.9 
we get the following result, which was obtained already by Berger [3]. 
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Theorem 5.1 (Robustness of EMD, [H Lemma 3]). Consider 

x = A(t)x, (5.6) 

with A <= (7(1, L(R' 1 )). Suppose the associated solution operator $ admits an EMD 
on I. T/ien i/iere exists 8 £ R>o swcft i/iai £/ie solution operator of any B G -B(A, <5) 
admits an EMD on I (wit/i i/ie same extremal projection as A). 

The last result is interesting from the following point of view. When imposed on 



R>o the inequalities (3.2) in Lemma 3.2 required for an EMD on R>o correspond to 



the definition of a so-called semistrong dichotomy of (5.6) on R>o (see |35|L The 
work by Vinograd [34, 35j yields that semistrong (exponential) dichotomies on R>o 
are robust only in the larger class of general exponential dichotomies, but not within 
semistrong dichotomies. That means that robustness of EMD can not be deduced 
from the classical asymptotic analysis but is a pure finite-time result. 



Analogously to the robustness investigation for linear processes in |subscction 3.3 



we now address the question of the stability radius for linear ordinary differential 
equations given on a compact time-interval. 



Definition 5.2 (Stability radius). Suppose Eq. (5.6) with A 6 C(I,L(R n )) generates 



an attracting solution operator on /. Then we define the stability radius of A by 

r(A) := inf {UBII^ ; B E (7(1, L(R n )), (A + B) is not attracting on 1} . 

To calculate the stability radius of a given A we make use of an elementary result 
which can be found in [10 and which specializes to the following result. 

Proposition 5.3 (cf. [10, Proposition 1, p. 2]). Consider 

x = A{t)x, 

with A £ C(I, L(R™)). Suppose the associated solution operator $ is attractive on I, 
i.e. $ admits an EMD on I with the trivial projection idjRn and A(R™, $) < 0. Then 
for any B G C(I,L(R")) with H-B^ < |A(R",$)| =: 8 the solution operator * of the 
perturbed ODE 

y = (A(t) + B(t))y 
satisfies A(R™, *) < A(R'\ $) + 8 = 0. 

In other words, — A(R n , $) > is a lower bound on the stability radius around 
A in (7(1, L(R n )). The fact that it is the stability radius follows directly from the 



well-known correspondence of shifted ODEs and weighted processes (cf. Eq. (3.4)), 
the definition of spectrum based on the weights / shifts and the fact that for operator 
norms induced by a vector norm the identity has norm 1. In summary, by a completely 
analogous argumentation as in the proof of |Thcorem 3.1l] wc find that 

y=(A(t)-\(R n ,<S>)id Rn )y 

is not attractive. Thus, we obtained that the stability radius of an attractive linear 
ODE given on I and the (pseudo-)stability radius of its associated solution operator 
coincide. 
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Theorem 5.4 (Finite- Time Stability Radius). Let A G C(I,£(R n )) and let the asso- 
ciated solution operator $ be attractive. Then 

r(A) = -A(R n ,$). 

Another consequence of EMD-robustness deals with linearizations. The next lemma 
states that the linearization depends continuously on the initial value of the trajectory 
along which we linearize. 



Lemma 5.5. Consider Eq. (5.1| where f G C^I x Z),R™) satisfies conditions for 
well-posedness. Let if denote the associated solution operator. Then the function 

D -> C(I,L(R")), X ^ dif(;<p(;,t-,x)), 

is continuous. 

Proof. Let t G I, x G D and e G R>o- First note that y ip(-,t-,y) G C(I, R n ) 
is uniformly continuous on any bounded subset U C D with x G U, i.e. there exists 
Si : R>o — ► R>o such that for any e G R>o and y G U we have 

||^(-, i) y)^ <e. 

By the C 1 assumption on / we have in particular that there exists some function 
62 ■ R>o ~~ ^ R>o such that for any e G R>o we have 

\<p(t,t-,x) - <p{t,t-,y)\ < 62(e) =► ||9i/(t,^,t_,a;)) < e. 

Combining the two continuity observations we obtain for any e G R>o and y G R n 

\X- V \ <6 1 (S 2 ( S ))^\\d 1 f(;<p(;t-,x))-d 1 f(;<p(;t-,y))\\ 00 < E. □ 

Consequently, robustness of EMD carries over to initial values. 



Corollary 5.6 (|H Theorem 5]). Consider Eq. ( |5.1| with f G C 1 (E x L», R") satisfying 
conditions of well-posedness. Let if be the associated solution operator and tp(-,t—,x) 
be a hyperbolic / attractive / repulsive reference solution. Then there exists 6 G R>o 
such that for any y G B(x,6) the trajectories <£>(•, i_,y) are, respectively, hyperbolic / 
attractive / repulsive. 



The next step is to show that Lemma 4.11 applies. Therefore, it remains to establish 
the necessary order of convergence for the linearization error. To this end, consider 
Eq. ([51} and assume that / G C°' 2 (L x R",R n ). By Taylor's Theorem, the estimate 
(5.3) on the nonlinear term g improves as follows: there exist s,5 G R>o such that 

sup J J^jH; xeB(0,6) \{0}, t Gil <£. (5.7) 
\x\ 
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Lemma 5.7. Let <p be the solution operator of (5.2) with f e C°' 2 (I x It™, It 71 ) 
and <fr be the solution operator of the linearization (5.4 1 along the reference solution 
ip(-,t-,x) = 0. Then the function 



fo, 

r\ i — ^ 



77 = 0, 



^sup{||<^(-,t_,t/) - $(-,t_)j/|| 0O ; yeB[0,77]}, otherwise, 
is continuous in 0. Moreover, it is of class 0(r] 2 ) for r\ — > 0. 



Proof. Let r/ g R >0 and y S It™, \y\ < rj. Integrating Eqs. (5.2) and (5.4), we calculate 
for any t El 



\p(t,t_,y)-$(t,t_)y\ = 



< 



+ 



+ 



A(s) (p(s,t_,y) - $(t,s)y) + g(s,cp(s,t_,y))ds 



\A(s)(<p(s,t-,y) - $(a,t-)»)l ds+ 



\g(s,(p(s,t-,y)} - g(s,$(s,t-)y)\ ds+ 



\g(s,$(s,t_)y)\ ds 



<m^\L+ L f) J \<p(s,t_,y)-$(s,t_)y\ds+ 
+ (<+ - t_)eC^rf, 



where C$ := ||$(- 



< oo and e > satisfies Eq. (5.7). With the abbreviation 



C{ri) := (t+ - <_)eC$7 7 2 e( 2 ll A ll-+ L /)( t +- t -) e 0( V 2 ), 
the uniform estimate from above and Gronwall's lemma we obtain 

B up{||v>(-,i-,iO - ^C-.*-)»lloo ! V e ^ ^fa) 

and the assertion is proved. 



□ 



Lemma 5.8. Let <p be the solution operator of Eq. (5.2) with f € C°' 2 (I x It™, It") 
and $ fee the solution operator of the linearization (5.4) along the reference solution 
ip(-,t-,x) = 0. Then the function 



fo, 



rj i ^ 



77 = 0, 

sup{||9oV(-,*-,2/)-^o*(-i*-)y|L; 2/ S -B[0,r?]}, otherwise, 



is continuous in 0. Moreover, it is of class 0(r] 2 ) for rj — > 0. 
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Proof. Let r\ € R >0 and y £ R", \y\ < ij. We calculate 

\do<p(t, t_, y) - do$(<, = |A(*)(yj(t, *_,»)- $(t, t_)lO + ff(t, ¥>(*, t_,y))\ 

< HAI^ | V (t, «_,!,) -$(t,t_)y| + 

+ |s(t^(t,«-,y))-s(t,*(t,t-)i/)| + 

+ | ff (t,$(M-)tf)| 

< Halloo C(^) + (MIL + £/)C(»7) + eC*^ 
<2P|| oo +L / )C(r 7 )+eC $?7 2 , 

where we used the notation from the proof of |Lcmma 5.7| This proves the assertion. □ 



In summary we have that the function m as defined in Eq. (4.4 1 is continuous in 
and hence |Theorem 4.12| |Thcorcm 4.13| and |Thcorem 4.15| apply. 

Remark 5.9. We want to comment on some issues concerning the application of the 
results in lscction 41 to ODEs. 

1. Concerning |Theorem 4.13| note that EMD-subspaces from the starting time- 
fibre {t-} x R™ evolve non- linearly under tp. Their extensions via tp considered 
as subsets in the extended state space I x R™ are C 1 -manifolds since ip(-,t-,-) 
is continuously differentiable (see, for instance, [31 Theorem 9.2]). 

2. Our version of |Thcorcm 4.13| applied to the ODE situation extends some previous 
work on that topic: in [24l [2TJ [5] the ODE is first extended to the real line and 
the desired manifolds are then obtained as the standard stable and unstable 
manifolds of trajectories that are hyperbolic on R in the sense that they admit 
an exponential dichotomy. Note that our proof does not restrict to Banach space 
norms induced by an inner product weighted with a symmetric, positive definite 
matrix, as it is done in [5] (cf. also Remark 4(h)], where it is stated that 
the presented proof does not yield EMD-like estimates with respect to arbitrary 
differentiable norms). In [15) an "intrinsic" proof is presented which works only 
in the case R™ = R 2 . There it is shown, that the domains of attraction and 
repulsion are not empty under appropriate conditions, but no assertions on the 
structure of contained subsets are proved (locally subspaces, manifolds, etc.). 

Since our hyperbolicity notion is based fundamentally on the monotonicity of the 
norm of trajectories, a finite-time conjugacy between the linear process $ and the 
general process ip should preserve the type of monotonicity of trajectories. Roughly 
speaking and supposing that the assumptions are satisfied for I, Theorem 4.16| can 



therefore be interpreted as a finite-time Hartman-Grobman-like theorem in the follow- 
ing informal sense: as demonstrated in (331 p. 546] the function 

H: I x R" -> R n , (t,x) i — ^ <p(t,t-,$(t_,t)x), 

with fibrewise inverse 

Hit,-)' 1 : R™ — s- R™, x h-> <!>(t,t_,ip(t_,t)x), 
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maps trajectories of <E> to trajectories of tp homeomorphically and is therefore a candi- 
date for a (nonautonomous) topological conjugacy between $ and ip with respect to the 
two zero reference solutions. Now consider restrictions of H to V a ($) and V u (&) (or 
compact subsets S C V s ($) and U C V u (&) considered as subsets of Gr(l,R n )). By 



Theorem 4.16 we obtain for y from the respective set with \y\ sufficiently small, that H 



preserves the monotonicity type of $(•,<_)?/, but, in general, not the exponential rate. 
It is unclear whether in general for hyperbolic processes there exists a neighborhood 
U C K™ and a (nonautonomous) topological conjugacy H (with respect to the two 
zero reference solutions) which preserves the monotonicity behavior of trajectories. By 
our method of proof, problems seem to arise towards the boundary of the stable and 
unstable cone, not to mention trajectories with initial values outside the two cones. 



6 Conclusions 

In this work we introduced an abstract framework for what is often called "finite-time 
dynamics" by adapting the notion of (invertible) processes to the situation of a compact 
time set I. Based on the notion of growth rates which we defined independently from 
any hyperbolicity notion we managed to give a coherent, unified and comprehensive 
presentation of the theory of linear analysis of finite-time processes on I. Due to 
the strong ambition to establish continuity of the involved functions, we obtained 
new and simpler arguments to prove partially known results, sometimes getting rid 
of unnecessary technical assumptions such as differentiability or the finiteness of the 
time set. 

Evidently and once more shown by our results, finite-time hyperbolic trajectories 
play an important role in the local analysis of finite-time dynamics. Therefore, to 
have sufficient, possibly computable criteria for hyperbolic solutions is necessary for 
the application of the hyperbolicity concept to real-life problems. For two-dimensional 
Hamiltonian systems Haller [53] proved the existence of a hyperbolic solution in the 
neighborhood of a curve of so-called instantaneous stagnation points under some ve- 
locity bound on the governing vector field, whereas Due & Siegmund proved recently 
the existence in the neighborhood of an approximate solution JTHj . In [T3] it is shown 
that row diagonal dominance of the linear right hand side implies an EMD on I. Other 
sufficient criteria are based on a partitioning approach [551 [JJ [T3] . 

We would like to emphasize that our way of reasoning is led by the idea to consider 
finite-time analysis as the analysis of "unextendible" systems, thereby restricting any 
argumentation to the time set I. This allows for a consequent use of continuity ar- 
guments. However, other methods of analysis and proof may become available and 
reasonable when approximating infinite-time dynamical systems on bounded time-sets. 



Acknowledgments 

The author is indebted to Maik Groger, Sascha Trostorff and Marcus Waurick for 
several enlightening discussions. He also would like to thank Stefan Siegmund for 



30 



introducing him to finite-time dynamics and for a couple of comments and Martin 
Rasmusscn for bringing [3] to his attention. 



References 

P.-A. Absil, R. Mahony, and R. Sepulchre. Riemannian Geometry of Grassmann 
Manifolds with a View on Algorithmic Computation. Acta Applicandae Mathe- 
maticae, 80(2):199-220, 2004. 

B. B. Aldridge, G. Haller, P. K. Sorger, and D. A. Lauffenburger. Direct Lyapunov 
exponent analysis enables parametric study of transient signalling governing cell 
behaviour. IEE Proceedings - Systems Biology, 153(6) :425-432, 2006. 

H. Amann. Ordinary differential equations: an introduction to nonlinear analysis, 
volume 13 of De Gruyter studies in mathematics. Walter de Gruyter, 1990. 

A. Berger. More on finite-time hyperbolicity. Boletin de la Sociedad Espanola de 
Matemdtica Aplicada, 51:25-32, 2010. 

A. Berger. On finite-time hyperbolicity. Communications on Pure and Applied 
Analysis, 10(3):963-981, 2011. 

A. Berger, T. S. Doan, and S. Siegmund. A Remark on Finite-Time Hyperbolicity. 
PAMM, 8(1):10917-10918, 2008. 

A. Berger, T. S. Doan, and S. Siegmund. A definition of spectrum for differential 
equations on finite time. Journal of Differential Equations, 246(3):1098-1118, 
2009. 

P. Bohl. Uber Diffcrentialungleichungen. Journal fur die Peine und Angewandte 
Mathematik, 144:284-318, 1913. 

F. Colonius, P. E. Kloeden, and M. Rasmussen. Morse spectrum for nonau- 
tonomous differential equations. Stochastics and Dynamics, 8(3) :35 1-363, 2008. 

W. A. Coppel. Dichotomies in Stability Theory, volume 629 of Lecture Notes in 
Mathematics. Springer- Verlag Berlin Heidelberg New York, 1978. 

C. M. Dafermos. An invariance principle for compact processes. Journal of Dif- 
ferential Equations, 9:239-252, 1971. erratum: ibid. 10 (1971), 179-180. 

Yu. L. Daletskij and M. G. Krejn. Stability of solutions of differential equations in 
Banach space. , volume 43 of Translations of Mathematical Monographs. American 
Mathematical Society (AMS), Providence, R.L, 1974. 

T. S. Doan, D. Karrasch, N. T. Yet, and S. Siegmund. A unified approach to 
finite-time hyperbolicity which extends finite-time Lyapunov exponents. Journal 
of Differential Equations, 252(10):5535-5554, 2012. 



31 



141 T. S. Doan, K. J. Palmer, and S. Siegmund. Transient Spectral Theory, Stable and 
Unstable Cones and Gershgorin's Theorem for Finite-Time Differential Equations. 
Journal of Differential Equations, 250(11):4177-4199, 2011. 

151 L. H. Due and S. Siegmund. Hyperbolicity and Invariant Manifolds for Planar 
Nonautonomous Systems on Finite Time Intervals. International Journal of Bi- 
furcation and Chaos, 18(3):641-674, 2008. 

161 L. H. Due and S. Siegmund. Existence of finite-time hyperbolic trajectories for pla- 
nar Hamiltonian flows. Journal of Dynamics and Differential Equations, 23:475- 
494, 2011. 

171 J- Ferrer, M. I. Garcia, and F. Puerta. Differentiable families of subspaces. Linear 
Algebra and its Applications, 199:229-252, 1994. 

181 P- Giesl and M. Rasmussen. Areas of attraction for nonautonomous differential 
equations on finite time intervals. Journal of Mathematical Analysis and Appli- 
cations, 390(l):27-46, 2012. doi: 10.1016/j.jmaa.2011. 12.051. 

191 I- C. Gohberg, P. Lancaster, and L. Rodman. Invariant subspaces of matrices 
with applications. Classics in applied mathematics. SIAM Society of Industrial 
and Applied Mathematics, 2006. 

201 G. Haller. Finding finite-time invariant manifolds in two-dimensional velocity 
fields. Chaos, 10(1):99-108, 2000. 

211 G. Haller. Distinguished material surfaces and coherent structures in three- 
dimensional fluid flows. Physica D, 149(4) :248-277, 2001. 

221 G. Haller. Lagrangian structures and the rate of strain in a partition of two- 
dimensional turbulence. Physics of Fluids, 13(ll):3365-3385, 2001. 

231 G. Haller. A variational theory of hyperbolic Lagrangian Coherent Structures. 
Physica D, 240(7) :574-598, 2011. 

241 G. Haller and A. C. Pojc. Finite time transport in aperiodic flows. Physica D, 
119(3):352-380, 1998. 

251 G. Haller and G. Yuan. Lagrangian coherent structures and mixing in two- 
dimensional turbulence. Physica D, 147(3-4) :352-370, 2000. 

26] T. Kato. Perturbation Theory for Linear Operators. Classics in Mathematics. 
Springer- Verlag Berlin Heidelberg New York, corrected printing of the second 
edition, 1995. Reprint of the 1980 edition, Originally published in Grundlehren 
der mathematischen Wissenschaften; Bd. 132. 

[27] T. Peacock and J. Dabiri. Introduction to Focus Issue: Lagrangian Coherent 
Structures. Chaos, 20(1):017501, 2010. 



32 



[28] M. Rasmusscn. Finite-time attractivity and bifurcation for nonautonomous dif- 
ferential equations. Differential Equations and Dynamical Systems, 18(l):57-78, 
2010. 

[29] K. Rateitschak and O. Wolkcnhauer. Thresholds in transient dynamics of signal 
transduction pathways. Journal of Theoretical Biology, 264(2) :334-346, 2010. 

[30] R. J. Sacker and G. R. Sell. A Spectral Theory for Linear Differential Systems. 
Journal of Differential Equations, 27(3):320-358, 1978. 

[31] S. C. Shadden, F. Lekicn, and J. E. Marsden. Definition and properties of 
Lagrangian coherent structures from finite-time Lyapunov exponents in two- 
dimensional aperiodic flows. Physica D, 212(3-4):271 -304, 2005. 

[32] S. Siegmund. Dichotomy Spectrum for Nonautonomous Differential Equations. 
Journal of Dynamics and Differential Equations, 14(l):243-258, 2002. 

[33] S. Siegmund. Normal Forms for Nonautonomous Differential Equations. Journal 
of Differential Equations, 178(2):541-573, 2002. 

[34] R. E. Vinograd. Exact bounds for exponential dichotomy roughness I. Strong 
dichotomy. Journal of Differential Equations, 71(1):63-71, 1988. 

[35] R. E. Vinograd. Exact bounds for exponential dichotomy roughness III. 
Semistrong dichotomy. Journal of Differential Equations, 91(2):245-267, 1991. 



33 



